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Abstract. 

Results from the Wilkinson Microwave Anisotropy Probe (WMAP), Atacama 
Cosmology Telescope (ACT) and recently from the South Pole Telescope (SPT) have 
indicated the possible existence of an extra radiation component in addition to the well 
known three neutrino species predicted by the Standard Model of particle physics. In 
this paper, we explore the possibility of the apparent extra dark radiation being linked 
directly to the physics of cold dark matter (CDM). In particular, we consider a generic 
scenario where dark radiation, as a result of an interaction, is produced directly by 
a fraction of the dark matter density effectively decaying into dark radiation. At an 
early epoch when the dark matter density is negligible, as an obvious consequence, 
the density of dark radiation is also very small. As the Universe approaches matter 
radiation equality, the dark matter density starts to dominate thereby increasing the 
content of dark radiation and changing the expansion rate of the Universe. As this 
increase in dark radiation content happens naturally after Big Bang Nucleosynthesis 
(BBN), it can relax the possible tension with lower values of radiation degrees of 
freedom measured from light element abundances compared to that of the CMB. 
We numerically confront this scenario with WMAP+ACT and WMAP+SPT data 
and derive an upper limit on the allowed fraction of dark matter decaying into dark 
radiation. 
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1. Introduction 

Measurements of fluctuations in the Cosmic Microwave Background (CMB) power 
spectra [1] and light element abundances from Big Bang Nucleosynthesis (BBN) [2] 
have been corner stones in the era of precision cosmology. Though standard ACDM 
cosmology is a very good fit to present data from these measurements, there are also 
tantalising hints of physics beyond standard ACDM. One such example is the possible 
presence of an extra dark radiation component during the epoch of decoupling. Recent 
analyses of CMB data point towards the possible existence of one or more radiation 
component (s) [3-7] with no standard electromagnetic and electroweak interactions other 
than those predicted by the standard model of particle physics. This extra radiation 
needs to be dark in the sense that the presence of an extra photon like component 
would not only spoil the success of BBN, but also generate a chemical potential for 
the photon - something which is constrained by CMB observations. The indication of 
excess radiation arises mainly through the precise observation of less power in the smaller 
scales of CMB anisotropy spectra. It has also been confirmed that these hints for extra 
radiation are indeed 'real', insofar as not being a statistical ambiguity from the choice 
of confidence interval [8]. Recently, the evidence was bolstered through the possible 
indication of ultra light sterile neutrino states in the neutrino oscillation experiments 
[9-13] and also through the reactor neutrino anomaly [14]. 

The nature of the dark radiation component is a topic of much debate. Current data 
allow the dark radiation component to be comprised by both sterile neutrinos as well 
as active neutrinos with a nonthermal distribution (see e.g. [15]). If dark radiation is 
comprised by massless sterile neutrinos, we expect them to behave as relativistic particles 
with effective sound speed Cg ff and viscosity parameter c^ is satisfying c^ s = cl is = 1/3$. 
Possible deviations from these values could indicate nonstandard interactions in the 
neutrino sector [15,17-20]. Luckily, measurements of CMB anisotropics can help in 
constraining these parameters [21] and most analyses are consistent with the c^ ff = 
c 2. s = 1/3 (see e.g. [7,22,23]) - although [24] reported on finding c 2 cS < 1/3. The 
bottom line is that cosmological data is sensitive to the details of dark radiation and 
can help in predicting its nature. 

The helium abundance Yp is very sensitive to the expansion rate of the Universe 
(and hence the amount of radiation present) at the time when T ~ MeV. However, the 
evidence for extra radiation from BBN data is somewhat ambiguous. In some analyses, 
it is reported that one can accommodate one extra dark radiation component [25], 
while in other works it is concluded that there is no need for extra radiation during 
BBN [26,27]. The former of these studies, for example, found JV;| BN = 2.4 ± 0.4 [26]. 
The main reason for the confusion is that N e g is highly sensitive to how the helium 
abundance is treated in the analysis [28]. In general, from many data analyses, it 
remains a possibility that the central value for the number of relativistic degrees of 
freedom allowed by CMB data is higher than that of BBN. Taking CMB data alone, the 

| Check definitions in [16]. 
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estimate is A^ MB = 5.3±1.3 from the Wilkinson Microwave Anisotropy Probe (WMAP) 
7-year with Atacama Cosmology Telescope (ACT) data [29]. Combined WMAP and 
South Pole Telescope (SPT) data give a slightly lower value of A C ^ MB = 3.85 ±0.62 [30]. 
The addition of baryon acoustic oscillations (BAO) data and the measurement of 
the Hubble parameter Hq improves these constraints somewhat. It is found that 
n cmb = 4 56 ± 0.75 for WMAP+ACT+BAO+# [29] and A e ^ MB = 3.86 ± 0.42 for 
WMAP+SPT+BAO+i/o [30]. Taken at face value, the latter two results suggest ~ 2a 
evidence for extra relativistic species. The Planck satellite will dramatically increase 
the precision of the inferred value of AN e g ~ 0.26 [31] and should be able to find a 
mismatch (if there is one) between N™ B and N®® N at the level of 4 — 5cr [3] . 

To explain the apparent radiation excess, one can, of course, just add a weakly 
interacting neutrino-like fermion by hand. However, then the question remains of ex- 
plaining the origin of such a particle. For a recent particle physics model explaining the 
radiation excess with three flavours of light right-handed neutrinos, though, see [32]. 
See also [33-35] for another explanation of excess radiation during decoupling through 
the Bose-Einstein condensation of a coherently oscillating dark matter axion. We would 
like to point out that though LSND [36] and MiNiBooNE [37] indicate the existence of 
one or more eV scale sterile neutrinos [38] , which are excellent candidates for the excess 
radiation hinted by small scale CMB data, it is very hard to reconcile two eV scale ster- 
ile neutrinos as dark radiation with the large scale structure and other measurements 
[39] unless sterile neutrinos have other interactions [40-44]. So, it is highly possible 
that the dark radiation may be a result of dark sector physics. For instance, if dark 
matter decays into dark radiation, that can explain the dark radiation excess and its 
effect could be found in the structure formation of the Universe. Note that some hidden 
sector models [45-47] motivated by other issues of particle physics and cosmology can 
also provide extra AN e g. 

If it is indeed the case that there is a change (increase) in the number of radiation 
degrees of freedom between the epoch of BBN and CMB, that will be an extremely 
interesting and surprising result. From a theoretical point of view, some new physics 
has to set in at a low energy scale (T ~ eV). Recently, there have been a few interesting 
works in this line of thought [48, 49], where AA^ 61 ^ ^ AA^ MB . From a particle physics 
view point, this indicates that a particle (beyond the frame work of the standard model) 
has to decay [50] into an extra dark radiation component in between the epoch of BBN 
and photon decoupling. 

In this paper we propose a very simple mechanism where one naturally generates 
an extra radiation component when the Universe approaches the era of matter radiation 
equality and decoupling. The basic idea is to allow dark matter to interact with and 
decay into dark radiation: As the Universe approaches matter radiation equality (MRE), 
the density of dark matter starts to dominate the universal energy budget. As a result 
of the interaction, the densities of dark matter and dark radiation are proportional, 
hence the density of dark radiation increases as we approach MRE. In this scenario, one 
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would naturally see an increase in the dark radiation component after BBN but before 
decoupling. One extra advantage of this scenario is the fact that the decay naturally 
reduces the amount of dark matter in galaxies and clusters. As a consequence it may 
help to alleviate [51, 52] the well known small scale structure issues in ACDM cosmology 
- the problems with cuspy cores and overproduction of satellite galaxies in numerical 
simulations of structure formation [53-55] . We leave the details of this as possible future 
work. 

In this paper we solve for the dark radiation density as a function of redshift 
numerically and show that we can obtain AN e g 4 1 as the Universe approaches 
the epoch of photon decoupling. We confront this generic scenario with the present 
cosmological data. We take a model-independent approach, where we use WMAP with 
either ACT or SPT data§ to constrain the fraction of the dark matter density which is 
allowed to be converted into dark radiation. We show that one can easily find a viable 
region in parameter space where AA^ MB can be greater than AA^ BN by of order unity. 

The plan of the paper is as follows: In section 2, we quantify the production of dark 
radiation and solve for the background solution. We present an analytical expression 
of AAefj and plot its dependence on scale factor and the coupling between dark matter 
and dark radiation. In section 3, we derive the cosmological perturbation equations for 
our scenario and, in section 4, we discuss our main numerical results from a COSMOMC 
analysis using various datasets. We demonstrate that observations are consistent with 
AA c flf ~ 1 around decoupling and AA c fj ~ around BBN. In section 5, we present a 
specific model in which the dark radiation production can be realised and constrain the 
model parameters. Finally we conclude in section 6. 

2. Interaction between dark radiation and dark matter 

2.1. Background evolution 

If dark radiation belongs to the dark sector along with dark matter, an interaction 
between the two could be possible. A general coupling (at the background level) can be 
described by the energy balance equations 

Pdm + 3-fTpDM = —Q j 

Pdark + 3H(1 + U>dark) Pdark = Q, (1) 

where pdm and Pdark are the dark matter /radiation energy densities and H = a/ a is 
the Hubble rate, where a is the scale factor and an overdot denotes the derivative with 
respect to conformal time r. For our case of dark matter being converted into dark 
radiation then u>d ar k = -Pdark/Pdark = 1/3. The rate of energy transfer is given by Q - a 
positive Q denotes the direction of energy transfer from dark matter to dark radiation. 
A non-zero Q means that dark matter no longer redshifts exactly as 1/a 3 and also that 

§ We do not use ACT and SPT data simultaneously. This is because the observational fields slightly 
overlap and hence require a more detailed analysis of the combined noise properties. We thank Mark 
Halpern for pointing this out. 
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dark radiation does not redshift as 1/a 4 . It is important to remind ourselves that we 
require a covariant form for the energy momentum transfer Q. 

Several papers over the recent years have studied different forms of the energy 
transfer rate Q in the context of interacting dark matter-dark energy [58-68] . We adopt 
the covariant form of energy momentum transfer 4-vector introduced in [68] 

Qdm = TpDM , (2) 

where the form of interaction rate V depends on the details of the particle physics of the 
decay process. Many forms of T has been studied in literature, we adopt a simple case 
where T = aH, where a is a constant and H is the Hubble rate. As discussed in [68], an 
implicit assumption behind this form of T is that the interaction rate varies with time 
but not with space, which explains the presence of H in the place of the interaction 
rate in Eq. 2. This form of T can arise from different models of dark matter decay. In 
section 5, we demonstrate a model of dark matter decaying into dark radiation and we 
show that the above mentioned form of F can be easily realised in nature. 

For this form of the coupling it is easy to solve the background energy density 
equations [68] 

Pdm = PDM,o a ^ 3+Q ^ ) 

Pdark = Pdark^- 3 ^" + ( £ ) PDM^" 3 (a~ 3w ^ ~ a'") . (3) 

With u>dark = 1/3 the equation for pdark can be collected into two terms 

Pdark = /3a" 4 + (j^) PDM,oa" (3+a) , (4) 

where /3 is a constant. The first part behaves like a standard radiation density and the 
second part behaves like a fluid with an equation of state a/3. In the case of a weak 
coupling between dark matter and dark radiation, we require that a is small, which in 
turn leads to (3 ~ 0. In the following we only keep the second term in Eq. 4. This is 
further justified by the fact that the fraction of Pdark, which redshifts like 1/a 4 , will be 
subdominant to the fraction which redshifts like dark matter due to the expansion of 
the Universe. 

With this assumption we obtain 

Pdark _^ Q? _ Ct 

Pdm 3w d ark - a 1 - a 
As we will see later from our numerical analysis, a < 1, so the ratio reduces to 
Pdark/PDM = ol- In fact, in section 6 we study a model of dark matter decay for a 
specific interaction and decay mechanism where our assumption is realised. This means 
that the pure radiation-like component (the term with coefficient 0) is indeed absent. 

2.2. Calculation of AN eS 

In the standard cosmological scenario, it is a standard practice to define AiV cff by 

Prad = 1 + o JV cff | TFT J Pj > ( 6 ) 
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where the radiation density p rac j is given as a sum of the energy density in photons 
p 7 = (it 2 /15)T^ and standard model neutrinos. In the standard model this predicts 
TV™ _ 3_ 046 with T u /T y = (4/11) 1 / 3 . Any departure from the standard scenario is 
parameterized as N cS = + AN e g. 

One implicit assumption in the above definition is that the dark matter dilutes as 
1/a 3 . But in our case neither dark matter nor dark radiation dilutes in the standard 
way. This means we cannot simply use the definition above - a more appropriate model 
independent method is needed to compare the expansion rate H to that of the standard 
model Hsm and attribute the difference to AN c g. In our model 

3H 2 M 2 Pl = p DM ,o/a 3+a + T ^p D M,o/a 3+0 + p rcs t (7) 

I — a 

= Pdm,o/ a a + Prest j 

1 — a 

where p rest stands for the normal radiation and dark energy components. We then 
compare to the standard Hsu (with non-zero AN e s) to obtain 

5 Ajv -(|) 4 ^-T^"™> 3+ °-^ 3 - < 8 > 

Using the above definition we find that AN e s depends on the decay constant a as well 
as scale factor a. We show AN e g as function of a in the top-left panel of Fig. 1 for 
a = 0.02 and 0.04. At the time of decoupling this produces a AN e g of order unity. 

As already discussed, the contribution to AN c g arises mainly due to the faster 
expansion rate at the time of decoupling and MRE in our model. It is interesting to 
note that at an early epoch deep in the radiation-dominated era, when the dark matter 
density is negligible, the deviation from the standard expansion rate is close to zero, 
as expected. At late times - after z ~ 100 - the expansion rate again approaches the 
standard model expansion rate, as demonstrated in Fig. 1. This is because a — > 1 and a 
is still small - in other words the interacting dark matter starts to behave in the same 
way as in standard ACDM. 



3. Perturbations of dark radiation 



To determine the perturbations to dark radiation and dark matter we need to consider 
the Boltzmann equation for its distribution. For simplicity we will consider a massive 
dark matter particle decaying into a pair of massless daughter particles. Following 
standard practice (e.g. [69]) we expand the distribution function for each species j in 
terms of a zero-order component f® and a perturbation tyj 

fjix^q^n^r) = ff(q j} r)[l + ^j(x\ q j} m, r)] , (9) 

which depends on position x l , magnitude of momentum qj, direction and conformal 
time r. The phase space of each species obeys the Boltzmann equation 

A/i = ^ + ^d^ + df ± dq l + df l dri ± = fdfj\ 
dr dr dx l dr dqj dr dni dr \dr J c 
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Figure 1. For each panel we show the best-fit vanilla 6-parameter model from 
WMAP+SPT (black), then models with the same parameters but one extra relativistic 
species (dotted-red), a lower dark matter density of JIdm = 0.085 (dashed-blue, as 
opposed to Odm = 0.112) and a decaying dark matter model with a = 0.02 (dot-dash 
green) and 0.04 (dot-dot-dash magenta). (Top-left) AN c g as a function of scale factor. 
(Top-right) Hubble rate compared to the standard model. (Bottom-left) Effective 
(total) equation of state. (Bottom-right) Ratio of the gravitational potential <I> for a 
Fourier mode with k = 0.02 Mpc" 1 compared to the standard model. Horizon entry 
for this mode is indicated by the dashed vertical line. 



where is the collision term, which depends on particle interactions. 

At zeroth-order the Boltzmann equation for the dark matter distribution function 
can then be written as [70, 71] 

/dm = -«^/dm. (11) 

under the assumption that Eq. 5 is fulfilled. Upon multiplying by the proper energy 
6j = ^Jq? + a 2 m 2 and integrating over all momenta one obtains the same continuity 
equation for dark matter as in Eq. 1. 

Working in the synchronous gauge, we can now work out the equations of motion 
for the perturbations to dark matter and its decay product. For the perturbations to 
dark matter we write out the Boltzmann equation and the perturbation to the energy 
density by following the machinery described in Ref. [69]. In the end, the equations 
of motion for the dark matter perturbations reduce to the case of stable dark matter 
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particles - the only difference being the different distribution function specified by 
Eq. 11. This result was also obtained in Refs. [71-73] 

Following Ref. [69], for the massless decay product we integrate out the q 
dependence of the distribution function and expand the angular component in terms 
of Legendre polynomials, 

F '* ft - T) s ! -fSSSm 55 § H),(2 ' + I)F " (£> TWk ' ft) • (12) 

where \i = k ■ h and P n {[f) are the Legendre polynomials of order n. 
The Boltzmann equation can then be worked out to give 

2 • 4 /■ \ 

-Fdark + ikfiF dSuTk = -~h-~ I/1 + 677) P2W + H(l-a) (5 DM - Fd^k) ,(13) 



3 3 

where the expression for the density perturbation 5 = 5p/p in terms of phase-space 
integrals, and the definition of the synchronous gauge metric perturbations h and 77 can 
be found in Ref. [69]. 

Eq. 13 can be translated into a hierarchy of perturbation equations of motion for 
the dark radiation by inserting the expansion in Eq. 12 and collecting terms. The final 
results are 

4 2 • 

<Wk = - -#dark - ~h - H(l - «)(5 d ark ~ S DM ) , (14) 
^dark = k 2 ^-<5dark ~ Odark^ — H(l — «)6 l dark , 

8 3 4-8 

-Fdark2 = 2(T dar k = — 9 S ~ ~kFdaA3 + ~ h + -Tj — H(l — a)F dar k2 , 

15 5 15 5 



P 



k 



darkZ 



^dark(i-l) — (I + Ij-^dark (l+l) 



21 -L- \ I" v'-j-; \~ 1 -/-—uT^j H(l - a)F dark i , l>3 

where <5 d ark = ^darko, ^dark = 3/c/4F dark i and cr d ark = ^dark2/2. This is an infinite 
hierarchy so we need to truncate the hierarchy at some Z max , for which we choose [69] 

Fu(l m ^ + l) ~ "T Fvlmex - -^(/max-1) • (15) 

These equations are equivalent to those in Ref. [71] for their choice of decay variables. 

We implemented these equations in a modified version of CAMB [74]. This amounts 
to: (1) Changing the scaling behaviour of dark matter to a~( 3+ct ); (2) Modifying the 
background evolution to include an additional component whose energy density scales 
like dark matter but with u>d ar k = 1/3 ; (3) Implementing the hierarchy of perturbation 
equations, which are similar (with the exception of the final terms in Eq. 14) to the 
existing massless neutrino perturbation equations. 



4. Results 



In Fig. 2 we show the temperature power spectrum each of the models listed in Fig. 1. 
The general features can be understood qualitatively by the following 
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Figure 2. Temperature power spectrum for the models listed in Fig. 1. 

• For models with an extra relativistic species, lower CDM and decaying dark matter 
(DDM) the acoustic peaks are shifted to smaller scales. The angular scale of the 
peaks is set by the ratio of the sound horizon at decoupling to the angular diameter 
distance to decoupling. For a flat Universe this is approximately the ratio of the 
conformal time at decoupling to that today, i.e. 9a ~ t^ cc /t . For extra relativistic 
species the increased Hubble rate at early times decreases Td cc , while tq remains 
similar. For lower CDM both Td ec and To increase, with the relative increase in To 
compared to the standard model greater. For DDM both Td ec and To decrease, with 
the relative decrease in Td ec greater. 

• The first two peaks are noticeably enhanced in models with an extra relativistic 
species and lower CDM, but are suppressed in the DDM model. This arises from the 
driving effect (modes entering the horizon during the radiation era are enhanced 
due to the decay of gravitational potentials) and (for the first peak) the early 
Integrated Sachs- Wolfe (ISW) effect. The total (effective) equation of state of the 
Universe for each model is shown in Fig. 1. For extra relativistic species and lower 
CDM, radiation domination (u> cff rs 1/3) is extended, while for the decaying dark 
matter model the Universe actually departs from radiation domination faster. This 
is somewhat opposite to what one might expect, since dark matter is decaying into 
dark radiation, but is a consequence of fixing f2 DM today to be the same as the 
standard model and it scaling as a~^ 3+a \ The result of this can be seen in Fig. 1, 
by plotting the gravitational potential $ for a mode with k = 0.02 Mpc -1 , which 
enters the horizon around a « 10~ 4 . Potential decay is suppressed in the DDM 
model by the time of decoupling. 

• The first peak is also affected by the early ISW - here the potential can still decay 
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after decoupling as the Universe is not completely matter dominated. For DDM the 
total equation of state is closer to the standard model at decoupling than in models 
with an extra relativistic species or lower CDM, so the early ISW contribution is 
smaller. There is, however, a late time ISW effect, resulting in more power on large 
scales. This is because the total equation of state w c s ^ 0, and also because the 
radiation decay product (whose energy density is q;Pdm) contributes a source of 
pressure and anisotropic stress. 

• On small scales there is reduced power in the damping tail of the CMB for both 
extra relativistic species and DDM. This is due to the increased expansion rate 
prior to decoupling, resulting in higher diffusion damping. The opposite occurs for 
lower CDM due to the decreased expansion rate. 

In order to confront the DDM model with observations we perform parameter 
estimation using a modified version of the COSMOMC package [75]. For our analysis 
we use data from the 7-year WMAP release [1], the 148 GHz 2008 ACT data [29] and the 
150 GHz 2008/2009 SPT data [30]. Due to the small overlapping sky coverage between 
ACT and SPT we consider WMAP + ACT and WMAP + SPT independently. We use 
software provided by each team to compute the likelihood of cosmological models. 




0.020.040.060.08 0.1 0.120.14 0.02 0.04 0.94 0.96 0.98 1 1.02 1.04 




Figure 3. Marginalized parameter constraints for WMAP + ACT (solid black) and 
WMAP + SPT (dashed red). 
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To parameterize these models we fit for 7 parameters, imposing the flatness 
condition i\ = 0||: the baryon density Q^,h 2 , cold dark matter density Odm/j 2 , Hubble 
parameter Hq = 100 h Mpc -1 kms" 1 , optical depth to reionization r, and the amplitude 
A s and spectral index n s of initial fluctuations. In addition we fit for the dark matter 
decay constant a, imposing a prior that a > 0. Since ACT and SPT observe much 
smaller scales than WMAP, marginalisation over foregrounds is also required, since 
these contribute to the small scale temperature power spectrum. We follow the same 
procedure as in the ACT and SPT analysis, marginalising over a combined thermal 
and kinetic Sunyaev-Zeldovich (SZ), and a clustered point source, template, together 
with a Poisson [Ci = const) point source component. The reader is referred to these 
references for more details on the templates used. This brings the total number of 
parameters fitted to 10 (7 cosmo logical and 3 secondary foreground parameters). We 
used the lensed theoretical CMB spectra from CAMB in our fits, since lensing is favoured 
by WMAP + ACT/SPT at the level of several a. 

The results of our analysis are shown in Fig. 3 and Fig. 4, where we show 
marginalised 1-dimensional and 2-dimensional parameter constraints. There is no 
preference for a non-zero decay constant, with a 2a upper-limit of a < 0.027 for WMAP 
+ ACT and a < 0.028 for WMAP + SPT. This likely stems from the different effect 
on the CMB spectrum for DDM than for extra relativistic species, in particular the 
suppression of the first peak instead of an enhancement. There is some degeneracy with 
other parameters, in particular the spectral index and Hubble parameter. Repeating 
our analysis for the standard model (a = 0) we find, for example, WMAP + SPT gives 
n s = 0.964±0.011 and h = 0.706±0.021, while allowing a to vary then n s = 0.976±0.015 
and h = 0.765 ± 0.052. 

5. Origin of interaction: A phenomenological model 

An empirical form of energy transfer has been introduced in most studies [58-68] where 
interactions between different cosmological sectors has been considered. Till now we 
also have chosen an empirical form r = a H for our model and used it to confront with 
data. For this form of coupling we have shown that the ratio of the energy densities of 
dark radiation and dark matter Pd&rk/pDM = a is practically constant in time for small 
a. Hence, putting an upper bound on a basically gives us an estimate of how large a 
fraction of dark matter is allowed to decay into dark radiation, obeying all cosmological 
constraints. 

Though the main goal of this paper is not a dark matter model, we will now discuss 
a specific phenomenological model where one can analytically derive the energy transfer 
equations between dark matter and dark radiation. We show that a time independent 

|| The effect of leaving S\ as a free parameter has been investigated in e.g. [56, 57]. Interestingly, in a 
model with two sterile neutrinos the cosmological constant seems to be ruled out at 95 % confidence 
level. Furthermore, models with sterile neutrinos seem to prefer w < — 1 for the dark energy equation 
of state. 
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Figure 4. Marginalized constraints on the 7 fitted cosmological parameters for WMAP 
+ ACT (solid black) and WMAP + SPT (dashed red). Likelihood contours show the 
68% and 95% confidence levels. 



a can indeed emerge in a phenomenological model. Note, however, that our numerical 
results are not limited to this specific model. Any dark matter model with a coupling to 
dark radiation and where the fraction of dark radiation to dark matter does not change 
much in the course of a Hubble time will be subject to the constraints from section 4. 

We consider a coherently oscillating scalar field which plays the role of CDM [76-78]. 
This has a Yukawa type coupling to a nearly massless dark fermion ipd (£ D A0VW^)- 
This type of CDM can in principle decay parametrically into dark radiation and the 
situation is very similar to the fermionic preheating scenario in the context of inflation 
[79]. However, the energy scale which we are considering here is much lower compared 
to that of inflation. We refer our readers to Refs. [76-78] for details of the dark matter 
decay process in this scenario. Here we present a brief review about the basic mechanism 
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of dark radiation production from CDM and finally we put constaints on the model 
parameters on the basis of the numerical results obtained in section 4. 

By adopting the results of fermionic preheating [77, 79] in an expanding background, 
the comoving number density of dark fermions can be found by solving the well known 
Mathieu equations 

X'l + [k 2 + (m + ^qf) 2 - iy/qf]X k = 0, (16) 

where the resonance parameter q = A^g/m^, <pof(t) is the background solution for the 
time evolution of the oscillating scalar field, k = k/m^ is the dimensionless fermion 
mass, and fh = m^/m^. These three parameters completely determine the parametric 
production of fermions. We consider the oscillation of the field with the usual quadratic 
potential V = |m 2 </> 2 , which is a good first order approximation around the minima 
of any potential. The term (fh + y/qf) can be thought of as an effective mass of the 
fermion. As the scalar field oscillates, the effective mass itself will oscillate around zero 
and the parametric production of fermions is enhanced when the effective mass crosses 
zero. It can be shown numerically that rik(t) oscillates and due to Pauli blocking its 
maximum value never crosses unity. 

The expansion of the Universe has not been taken into account in the discussion 
above. As the mass of the oscillating scalar is very low, for a complete treatment, 
however, we have to consider the expansion of the Universe. When the expansion is 
taken into account, the resonance parameter q = A 2 4>(t) 2 / 'm^ becomes time-dependent 
and the periodic modulation of the comoving number density does not hold any more. 
It has been shown [77] (though in a different context of neutrino cosmology) that 
the parametric production of dark radiation happens as long as the time dependent 
resonance parameter satisfies q(z) ^> 1. In this regime, the produced dark radiation 
density takes a very simple form pdark = 87rA 2 pDM5 which gives a in terms of model 
parameters a = 8tt\ 2 . The only requirement for the above relation is q ^> 1, which 
translates into 

2*^(1 + *) 3 »1, (17) 

which we assumes holds true from BBN till the present epoch denoted by the '0' 
superscript. 

So far the dark matter mass has not been constrained in our analysis. This is 
because the dark matter mass does not enter explicitly into the perturbation analysis 
and thus our numerical results do not constrain the dark matter mass directly. One 
can, however, place an upper bound on the dark matter mass for this specific model 
of parametric decay. This is possible from the requirement > H - otherwise the 
Hubble friction would prevent the field from oscillating coherently and will not allow it 
to behave as CDM. From the constraints on CDM matter power spectra we know that 
dark matter has to be present in the Universe at least couple of e-foldings before MRE 
- otherwise there would be too much suppression in the linear matter power-spectra 
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on small scales [76]. So to get an estimate we assume the scalar started to oscillate 
coherently when the temperature of the Universe was around T = T osc ~ 100 eV. This 
choice keeps us out of the conflict with constraints from linear matter power spectra 
measurements from Lyman-a and SDSS data [80,81]. This, in turn, constrains the dark 
matter mass ^> H(T osc ). In Fig. 5 we show the allowed region in the (m^, A) plane 
which satisfy all of the above three constraints namely: a) To satisfy the parametric 
production of dark radiation; b) To be consistent with our numerical results (upper 
bound on a) from section 4; c) To obey the condition for the coherent oscillation setting 
in before MRE. 




0.00 0.02 0.04 0.06 0.08 0.10 
A 



Figure 5. Allowed region in the (ro^, A) plane for a coherently oscillating scalar dark 
matter decaying into dark radiation. The blue area is excluded from the requirement 
that the resonance parameter is not high enough to produce dark radiation through 
parametric resonance. The purple region is excluded from the upper limit on the 
fraction of dark matter decaying into dark radiation taken from our numerical result 
in section 4. The lower grey area is excluded from the requirement 3> H(T osc ). 



6. Discussion and conclusion 

The evidence for the existence of dark radiation at the CMB epoch is intriguing. If 
future experiments find a mismatch between the radiation content of the Universe at 
the epoch of BBN and decoupling, the production of dark radiation may be a late-time 
phenomenon in the cosmic history which took place some e-foldings after BBN. Future 
surveys like Planck will measure the effective number of radiation degrees of freedom 
with an accuracy of AN e ff = 0.026 [31] and will also be able to probe if extra radiation 
has been produced after BBN at all, so that AN^ N ^ AA^ MB . 

In this paper, we have shown that if dark radiation is produced from dark matter 
decays, the Universe naturally gets populated with an extra radiation component after 
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BBN but before photon decoupling. The reason is that as the Universe cools the dark 
matter density increases and, as a result, so does the dark radiation produced from 
it. We have constrained the fraction of dark matter which is allowed to be converted 
into dark radiation using the WMAP7 + ACT and WMAP + SPT data. We find 
an upper bound on this fraction using a COSMOMC analysis and show that it is 
possible to get an increase in N cS by of order unity as the Universe approaches the 
epoch of photon decoupling. However, the effect on the temperature power spectrum is 
somewhat different than adding in extra relativistic species by hand, most noticeably 
in the suppression of the first acoustic peak. For this reason, if Planck confirms the 
mismatch between Ai\Qjj B and AN^P^ , it remains to be seen how well the decaying 
dark matter model fits data. 

As a phenomenological example we have presented a model of dark matter decay 
and calculated the decay rate as a function of the coupling between CDM and dark 
radiation. 

Dark matter decaying into dark radiation could also have important implications 
for cosmological structure formation. Very recently the observations of high redshift 
massive galaxy clusters [82] has put ACDM cosmology under stringent constraints [83- 
85]. In fact, the presence of extra dark radiation during the CMB epoch may play a role 
in resolving these issues [86]. Note in this context that in Fig. 1 we have showed that 
the decay of the gravitational potential during decoupling is suppressed in the decaying 
dark matter model. This may boost early structure formation as discussed above. The 
detailed study of this effect in the context of our model is beyond the scope of this paper 
and we leave it for future work. 
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